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RECONSTRUCTING JACOBI MATRICES FROM THREE 

SPECTRA 

JOHANNA MICHOR AND GERALD TESCHL 


Abstract. Cut a Jacobi matrix into two pieces by removing the n-th column 
and n-th row. We give necessary and sufficient conditions for the spectra of the 
original matrix plus the spectra of the two submatrices to uniquely determine 
the original matrix. Our result contains Hochstadt’s theorem as a special case. 


1. Introduction 


The topic of this paper is inverse spectral theory for Jacobi matrices, that is, 
matrices of the form 


(i.i) 
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i)i «i \ 

a\ b 2 a,2 
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on -i 


on -i 

bN J 


This is an old problem closely related to the moment problem (see (b| and the 
references therein), which has attracted considerable interest recently (see, e.g., [[j] 
and the references therein, (2) , ||, ||). In this note we want to investigate the 
following question: Remove the n-th row and the ?r-th column from H and denote 
the resulting submatrices by H (from b\ to b n - 1 ) respectively H + (from b n+ 1 to 
6 jv)- When do the spectra of these three matrices determine the original matrix HI 
We will show that this is the case if and only if R_ and H + have no eigenvalues in 
common. 

From a physical point of view such a model describes a chain of N particles 
coupled via springs and fixed at both end points (see [||, Section 1.5). Determining 
the eigenfrequencies of this system and the one obtained by keeping one particle 
fixed, one can uniquely reconstruct the masses and spring constants. Moreover, 
these results can be applied to completely integrable systems, in particular the 
Toda lattice (see e.g., ||). 


2. Main result 

To set the stage let us introduce some further notation. We denote the spectra 
of the matrices introduced in the previous section by 

(2.1) a(H) = {X j }f =1 , = a(H + ) = 
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Moreover, we denote by (pj)^ 1 the ordered eigenvalues of H~ and H + (listing 
common eigenvalues twice) and recall the well-known formula 

, , -i 

(2.2) g(z, n) = -- - 


UfMz ~ A i) z - b n + a n m +(z, n) + a^ l _ 1 m-(z, n)' 


where g(z , n) are the diagonal entries of the resolvent (H — z)~ l and m±(z , n ) are 
the Weyl m-functions corresponding to iJ_ and H + . The Weyl functions m±(z , n) 
are Herglotz and hence have a representation of the following form 


(2.3) 
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With this notation our main result reads as follows 
Theorem 2.1. To each Jacobi matrix H we can associate spectral data 

(2-5) {Ai}f =1 , (/^)f = -\ 

where (jj = +1 if pj G <j(H + )\cr(H_), <jj = —1 if pj G a(H_)\cr(H + ), and 


( 2 . 6 ) 


a n a t ~ 




a i a i+ a i-i a k 


*fl*j = Hk =t4- 

Then these spectral data satisfy 

(i) Ai < pi < X 2 < P 2 < • • • < Ajv, 

(ii) aj = a j+ 1 G (-1, 1) if Pj = p j+1 and aj G {±1} if Pj ± p t for i ± j 

and uniquely determine H. Conversely, for every given set of spectral data satisfying 
(■ i ) and (ii), there is a corresponding Jacobi matrix H. 

Proof. We first consider the case where H- and H + have no eigenvalues in com¬ 
mon. The interlacing property (i) is equivalent to the Herglotz property of g(z,n). 
Furthermore, the residues a~ can be computed from (EM) 


N—r 


rij=i ( 2 A?) ,2 

- 1 - - -a?-— = z-b n -a„ > 

n n—lf —\ y-rN— n/ 4-\ 11 n / v 
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(2.7) 
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and are given by a i = a n f 1 f3 i , where 


( 2 . 8 ) 


0T = - 


n;Li(M--A,) 


n— 1 


Similarly, af = a„ /3 ; , where 


11^7 - urmiTte - tf) 
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Hence m±(z, n ) are uniquely determined and thus H± by standard results from the 
moment problem. The only remaining coefficient b n follows from the well-known 
trace formula 


N 


n— 1 


N—n 


( 2 . 10 ) 


b n = tr (H) - tr(H_) - tr(tf+) = A i ~ ~ 




j =i 


fc=l 


z=i 


Conversely, suppose we have the spectral data given. Then we can define a n , 
a n - 1 , b n , a~p, oif as above. By (i), aj. and al are positive and hence give rise 
to H±. Together with a n , a n _i, b n we have thus defined a Jacobi matrix H. By 


construction, the eigenvalues /xjT, /r ; + are the right ones and also (2^) holds for H. 
Thus A j are the eigenvalues of H } since they are the poles of g(z, n). 

Next we come to the general case where /x JO = = (if (= A ;o ) at least for one 

jo- Now some factors in the left hand side of (2.7) will cancel and we can no longer 
compute /3^ o , 0^, but only 7 j 0 = 00 o + /?£. However, by definition of Oj 0 we have 


( 2 . 11 ) 


_ 1 ~ <Jjo 
Pk 0 o 


I301 


AI = 


1 


■ 730 - 


Now we can proceed as before to see that H is uniquely determined by the spectral 
data. 

Conversely, we can also construct a matrix H from given spectral data, but it is 
no longer clear that A j is an eigenvalue of H unless it is a pole of g(z, n). However, 
in the case A J0 = /j,^ = nf we can glue the eigenvectors of if_ and H + to give an 
eigenvector corresponding to A J0 of H. □ 


The special case where we remove the first row and the first column (in which 
case H_ is not present) corresponds to Hochstadt’s theorem Q. Similar results for 
(quasi-)periodic Jacobi operators can be found in Q. 
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